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Abstract 

In realistic model-building, spinor fields with various masses are present. During inflation, spinor 
field may induce gravitational waves as a second order effect. In this paper, we calculate the 
contribution of single massive spinor field to the power spectrum of primordial gravitational wave 
by using retarded Green propagator. We find that the correction is scale-invariant and of order 
H 4 /Mp for arbitrary spinor mass Additionally, we also observe that when i7u > H , the 

dependence of correction on m^/H is nontrivial. 
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I. INTRODUCTION 


General relativity strongly suggests the existence of gravitational waves. Through the 
study of gravitational waves, we might get rich information about new physics in the very 
early universe directly since there may be footprints left in it if the corresponding processes 
happened energetically. Recently, besides B-rnode polarization of CMB jl |, [2f|, gravitational 
waves searching are proposed or underway on smaller scale by interferometers like VIRGO, 
DECIGO and LISA, which make the studies relevant with gravitational waves acquire in¬ 
creasing attentions. 

The primordial gravitational wave may be produced from quantum fluctuations of vac¬ 
uum during inflation js| , ^] , { 5 ] , which is scale-invariant and is regarded as a smoking gun for 
inflation. The primordial gravitational waves may record some physical information during 
inflation, which can hardly be encoded by scalar perturbation, e.g. a step-like variation of 


primordial gravitational waves speed will result in the oscillating tensor spectrum 


6]. Thus, 


every possible way of affecting gravitational waves spectrum is significant and worth study¬ 


ing. There are also ot 


e.g. cosmic defeats 


fl 


rer mechanisms proc 
8|, bubble collisions 


ucing gravitational waves in the early universe, 


1Q|. 


Back in the very early stage of universe, there must be numerous fields besides the inflaton 
he scalar fields other than inflaton can only generate isocurvature 


field. In linear order, 


perturbations, se e |lllj . 12] for the effect of massive isocurvaton on the density perturbations, 
see also (l3].fl^|.jl5|. When the action is expanded to third order in perturbations, one may 
see that the scalar fields can source gravity at quadratic order. This motivates the mechanism 


of inducing gravitational waves by introducing a spectator scalar field 


16], see also 171, 18]. 


Spinor fields with various masses should also be present in early universe, as a conse¬ 
quence of realistic model-building. During inflation, these spinor fields look like spectators, 
i.e. they do not affect the inflationary dynamics at neither the background nor the linear 
order fluctuations level. Thus, similar to scalar spectators, they may also source gravity at 
quadratic order. In addition, there are also some studies about gravitational waves sourced 
by production of particles see e.g. 19], ( 20 I , by production of fermion see { 21 , 221. 


There are mainly two kinds of methods of calcu 
trum. One is loop correction with in-in formalism 


ating the corrections to the power spec- 


231]. The 


to the scalar power spectrum was studied comprehensively in 


oop correction of a Dirac field 


24]. In Ref. 25], we have calcu- 
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lated the loop correction of Dirac field to the gravitational wave power spectrum. The other 
is inducing gravitational waves as a second order effect, using retarded Green propagator, 
and regarding the spectator fields as sources. 

In this paper, we will calculate the contribution of massive spinor field during inflation 
to the gravitational wave power spectrum by using retarded Green propagator. We find 
that the correction is scale-invariant and of order H^/Mp for arbitrary mass m^. When 
m i> > H, the dependence of correction on m^/H is nontrivial. In Section II we give a 
brief introduction to the background dynamics of spinor held in inhationary cosmology. In 
Section III, we calculate the primordial gravitational wave spectrum induced by spinor held 
and give our result. The conclusion is given in Section IV. In Appendix A, we show the 
calculation process of spin sum for massive fermions, which is necessary for our computation 
of the massive spinor held’s contribution to the two point correlation function of tensor. 


II. DYNAMICS OF DIRAC FIELD 


In Ref. 


25| . we have described in detail the dynamics of a spinor held which is 


coupled to Einstein gravity, see also Refs. 


26 


27] for detailed introduction and 


minimally 


for the 


cosmology with spinor helds. In this section, we will briehy review it. 

Taking into account the spatial dependence of the Dirac held, the equations of motion 
for it is given by 


3 d ’ 

i7°(V> + - Hi/j) - = 0 , 


( 1 ) 


where H is the Hubble parameter during inhation. As usual, the Fourier transformation of 
V>(x, t) is 




d 3 P e ip x [«p,s^p,sO) + P-p,sW-p, s (t) , 

S 


( 2 ) 


where s = ±1/2 is the spin of the Dirac held, and also the annihilation operators of the 
Dirac held satishes the anti-commutative relation. 

The Fourier modes X ps {t), W PtS (t) satisfy the equations of motion shown in Eq.([T]). 
Solving the equation of motion in comoving time coordinate r = f dt/a, we obtain 


X k ,±(r) 


W-TikT , (1) 

2(27ra) 3 / 2 x ^ ir 


(-hr) , 


( 3 ) 
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where is the first kind of the Hankel function and r = m^/H, and also we have made 
use of Bunch-Davies vacuum as initial condition to determine the coefficient of the Hankel 
function. 


III. INDUCED GRAVITATIONAL WAVE 


In Einstein gravity, the quadratic action of tensor perturbation and the interaction be¬ 
tween it and the Dirac field are 


s (2 > = [ (fcg - OM 2 ). 


(4) 


S (3 l = / cfix—T^hij 

ipt/ih / O 3 

r n 2 

= i / d 4 x — (V’7*(<9- 7 , i/’) + hij, (5) 

respectively. Then the equation of motion for h^, sourced by Dirac field at second-order, 
can be expressed in conformal time as 


h"- + 2—ti- • - d 2 h■ = ^f-T TT 


( 6 ) 


where Tj 7 = U 1J lm Ti rn , and Ty is dehned in Eq.(j5]), and the TT-projection tensor n ij lm = 
P™Pj — \PijP lm , while Pij is the projection operator dehned by 


P.. — A.. 
1 %3 U 13 


kikj 

Tkp 


(7) 


where fcj/|k| denotes the direction of the propagation of a plane wave. The solution of Eq. 
takes the following form 

hij(k, t) = J dr , G fc (r,r , )a(r / )Tj T (k,r'), 

and Gfc(r, r') is the retarded propagator solving the homogeneous transform of Eq. d 
Within a de Sitter background a(r) = Gk(r,r') reads 


( 8 ) 


Gk(r, r') = ——-[(1 + fc 2 rr') sin k{r — t') + k(r' — r) cos k{r — r , )]0(r — r'). (9) 

k 6 r 2 

We will use Eqs. (JU [9]) in the following process. 
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Using Eq.(j 8 ]) and plugging in the Fourier transformation of the Dirac held, we get the 
tensor spectrum as 

{hij(k,T)hij(k', t)) 

= -—^4 J dr'a(r')Gk{T , r') j dr"a(r")G fc (r, r") J d 3 pd 3 p'U ij i rn (k)U ijno (k') 

<$(p, T') 1 l d m ^{k - p, r') - d l p(p, t') 7 m ^(k - p, r') + (lo m)] 

$(p', r")7 n ^(k' - p', r") - <9>(p', r'V^k' - p', r") + (n O o)]). (10) 


To proceed we calculate 


a 3 (r / )a 3 (r // ) y^tr[7 m Xp/ ]S /(T / )Xp/ ; y(r // )7 n VUp )g (r // )VUp iS (r / )] 

-pr"y 2ir , r ' 4 (h) 

~ p ' p " 

|r(/i)| 4 


r 4 f77') / - v 'r'\- 2ir 

= 2(p' m p n + p' n p m - p' ■ pS mn ){I — 1 


+2 6 r 


m 


(2jr) s V 2 

2nrf—2irff2ir , /2ivr r2ir ff—2ir\ 

-(e T T +e T T J, 


where /i = 1/2 — ir, see Appendix [X] for details. Then we have 


( 2 *) 


/ / \ 2ir 
-pr 


-pr 


2 ir 


(ii) 


(hij{k, r)hij(k', r)) 
_/Gk(T, t') 
a 2 (r') 


M 4 


dr' ~ K /y/ J I dr" 

,P 3 r 1 J.O 1 


Gfc(r,r") 

a 2 (r") 


| d 3 p<5 3 (k - k') 


{—16A—[- - (p • /c ) 2 + -(p • A;) 4 ] - 8 Ap • p'p'{ 1 - (p • fc) ] + 8 fip 2 [l - (p • fc) L ]},( 12 ) 
P z 


where A = 


r 4 A) UpV 

(2nY 


2 

—2ir 


—pr 


—2ir 


+ 


rV) (-p'r'\ 2ir (- P r"' 2ir 


and i? = 


IUaOI ! 2trr^i—2ir ^n2ir 


T 


2nr^/2ir _//—2zr N 


■^j-rrzir _j_ r^-rzir^. 


(2?r) 8 V 1 - 

with 0 and defining x = , 


(2tt) 8 \ 2 j \ 2 
After marking the angle between k and p 


= = t, y = T,we have 


x 


dcosO = — dx, 

y 


(13) 


since cos 8 = - + ^,p~ P = 1 + ' 9 y " :r , where p is regarded as constant in the integral before we 
do the y- integral. When cos# increases monotonically from -1 to 1, or 6 runs backward from 
7 T to 0, x = \Jl + y 2 — 2y cos 6 decreases monotonically from 1 + y to |1 — y\. Exchanging 
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the 9 and p integrals, we get 


(h ii (k,r)h ii (k , ,r)) 
5 3 (k — k')2nk 5 f°° 


M 4 


fi+ y 

dy / dxxy 

’ J\i-y\ 


2 /) , 1 4 m o 1 — — 2/ 2 2 

—16—-cos 0 + - cos 9) — 8 - 1 — cos i 

a+2 2 ; 2 xy K 


■T 4 (/r) f pp‘ 

LW) 


x 


/ \ — 2zr 


rfr' ^^V - 2 *] 2 + 


r 4 (/i) / pp 1 


2ir 


/G fc (r, r') 


a 2 (r') 


(2vr) 


£21 ] [ / d T ' ~ K y/; . V 2 "] 2 


a 2 (r') 


+8(1 — cos 2 0) 


\mr- 

(2t r) 8 


—2-7rr 


+e 27rr [ dT , Gk [ T / T '\ ' 2ir 


a 2 (r') 


j dT '^ T '- 2,r 
' dT " G -^y T "- 2 \ 


'^„ G k(T, T") M 

dT ^v+ T 


Note that for r = 0, 


Gfc (0, r') = [- sin hr’ + fcr' cos fcr']. 


(14) 


(15) 


Now, we take both time integrals to cover the range — oo —y 0. To do the time integrals, we 
combine k and t', t" to form new variables z' = kr r , z" = kr" . Then we have 


<Mk)Mk')> 

<5 3 (k — k')27rif 4 f°° 


Mpk 3 


fl+y 

dy / dxxy 
Jd-y\ 


x 


—16—(- - cos 2 9 + - cos 4 9) - 8-———^(1 - cos 2 0)] 
x y 2 2 ’ 2xy K 

(Xp-™ { [ dz ' z ^r { _ sm z , + , W)f + ac 


L(2tt) 8 V 4 


+8(1 — cos 2 9) 


(2tt) 8 


—2'Kr 


dz'z' 2ir (— sin z + z cos z '\) / dz"z" 2ir (— sin z" + z" cos z "]) 


+e 


2'Kr 


dz'z /2ir (— sin z' + z' cos z'}) / dz” z” ir (— sin z” + z” cos z”]) 


(16) 


where we have used a(r)H = —Integrating over time and momentum, we get 


(Mk)Mk')) 

5 3 (k - k ’)H A r 2 8ir 0F(2 i + r)T(2 - 2ir)T 4 (± + ir)T 2 (-27r) 

-Re J 


(2ir) 7 Mpk 3 


(i + 2r) 2 r(l — 2 ir) 


[4r(i + r) cosh(7rr)] 2 Z|, 
(17) 


where 


Z = 


i + r Exp{—27rr[189 + 3r(—187* + 6r(—41 + 22 ir + 4r 2 ))]} 


—5 + 2 ir 


(—9 + 2ir)(3i + 2r)(5i + 2r)(7i + 2 r) 


(18) 
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During the integration we encountered a lot of incomplete beta function B{ A; a, b ) and the 
hypergeometric function 2 Fi(a,b; c; A), where A is a physical cut-off, and 


B{z-, a, b) = z a [- + ^ + 0(^ 2 , z 3 ...)}, 

a 1 + a 


abz 


2 Fi(a, b;c;z) — 1 + AA1 _|_ 0(z 2 , z 3 ...). 


(19) 

( 20 ) 


c 

We may omit those terms ~ A n ,n > 1 which should be renormalized away. Expect these, 
we do not require any approximation relevant with r. Thus Eq. (1170 is an analytical result 
exactly for arbitrary r, i.e. m^/H. 

The result involves some Gamma functions written in term of r. However, we can plot its 
r-dependence out easily. In Fig. (JT]) , we plot the Re{...} part of Eq. flTTlh From Fig. ([Tj) , we see 


</iij/iij > (2 tt) 7 Mp k 3 


H 4 


H 4 



</iij/iij > (2 w) 7 Mj. k 3 


H 4 




(b) 


< /ijj Itij > (2 tt ) 7 Mp k 3 


H 4 



FIG. 1: The r-dependence of Eg. (1171) . We see that when r > 1, the Dirac field may contribute 


negatively or positively to the power spectrum of gravitational wave, depending on the value of r, 
i.e. m^p/H. 


that the r-dependent part is < 1, so the contribution of Dirac field to the gravitational wave 
power spectrum is ~ F 4 /Mp, which is scale-invariant. During inflation, the power spectrum 
of gravitational wave is ~ H 2 /Mp. Thus the contribution of spinor field is suppressed by 
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the factor H 2 /Mp. Additionally, from Fig. (JT]), we also observe that when m^ > H, the 
r-dependent part has a nontrivial behavior, i.e. the Dirac field with different mass may 
contribute negatively or positively to the gravitational power spectrum which manifests 
itself from the negative or positive fermion loop correction. 

The small mass limit r —>■ 0 of Eq. ffTTli is interesting, since this case provides a maximal 
correction, which is 


(hij (k) hij (k )) | r _ 


5 3 (k - k ')H 4 f 8tt 5 
(2vr ) 7 M 4 k 3 V '~75 


( 21 ) 


For a check, we follow Ref. 29], regulate the momentum integral by momentum cut-off, 
and replace the momentum integral and time integral as 


d 3 k —> 


C Aa(£l) 


d 3 k , 


r-r"(l+f) 


T ->• 


dr" 


dr ', 


( 22 ) 


where t' is the earlier one in time and A is a physical cut-off. Since the associative cut-off 
in comoving momentum is time dependent, therefore we first integrate over momentum and 
then integrate over time. After taking the small mass limit r —> 0, we have 


( kjj (k) hjj (k )) r— >-0 


<5 3 (k — k')A 5 
480 tt 4 A; 3 FM|(1 + f ) 4 

<5 3 (k - k!)H 4 /64tt 2 (56 + 44f + 18# + 3#) 

(27r) 7 M|A 3 V 75(l + f) 4 (2 + f) 3 

64tt 2 A 3 (4 + 28f + 78# + 56# + 44# + 18# + 3#) 
75 F 3 (l + f) 4 (2 + f) s 

647t 3 A 647t 3 A 3 87t 3 A 5 \ 

+ 63 H ~ 105 H 3 + 225F 5 (1 + f) 4 /’ 


where the first line on the right-hand side of Eq. (j23l) comes from the R-part of Eq. (fT2lh 
while the remaining parts of Eq. (fTTfll come from the A-part of Eq. (fEZD . We now leave 
behind the divergent terms ~ A n ,n > 1 assuming that they can be reabsorbed by counter 
terms in renormalization process, and we get 


S 3 (k - k')H 4 64tt 2 (56 + 44f + 18# + 3#' 


(hij (k) hij (k )) r _ 


(27 r) 7 M 4 k 3 75(1 + f ] 

<5 3 (k — k')i7 4 64 x 56 tt 2 


l (2 + 


(2tt) 7 M 4 A; 3 75 


(24) 


The above result shows the contribution induced by the Dirac field, in the small mass limit, 
is of the same order with Eq. (l2Tjh 















Until now, the impact of quantum field renormalization on cosmological correlation func¬ 
tions is controversial. Some authors argued that the influence of renormalization can signif¬ 
icantly change the predictions of slow-roll inflation for both scalar and tensor power spectra 


30h32|. While some authors debated that the standard expression is valid and the UV reg- 


ularisation leads to no difference in the bare power spectrum |33 h 371] . Though we focused on 


the effect of spinor held on the tensor power spectrum, similar case may also occur. However, 
the effect of spinor held is next-to-leading order, so we expect that available renormalization 
methods in literature do not alter our result in order of magnitude. 


IV. CONCLUSION 


Gravity, being the weakest interaction, might provide abundant information about new 
physics in the very early universe directly. Due to the weakness in strength of gravity, it is 
very hard to detect the signals of gravitational waves. However, forthcoming experiments 
might soon be able to catch some sight of it. This motivates us to explore possible sources 
of primordial gravitational waves. 

In this work, we studied the contribution of massive spinor held during inflation to the 
gravitational wave power spectrum by using retarded Green propagator. As an important 
step, we show the calculation process of spin sum for massive fermions in APPENDIX A. 
We worked out analytically the result for arbitrary mass m y, and found that the correction 
is scale-invariant and of order H^/Mp. Additionally, we also observe that when the mass 
of spinor held m y > H, depending on the value of my, the spinor held may contribute 
negatively or positively to the gravitational power spectrum. 

The result obtained is actually far below the observational level, which is reasonable since 
they are of the second order in perturbation. And of course, we are only limited to the case 
with single spinor held, if there are lots of spinor helds, our result will multiply. In addition, 
it is also interesting to consider non- minim a,! coupling to gravity, e.g. non-minimal derivative 
coiipling curvaton [3a, |39j, and things might be different. We leave this task to future work. 
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Appendix A: Trace 


Equation Eq. (TT0|) shows that the calculation of two point correlation function needs to 
sum over spins of massive fermions at different times. However, massive fermions in curved 
space are no longer conformally flat, and we can not use the spin sum formula from flat 
space. Here, we show the calculation process of this spin sum for massive fermions in great 
details. Despite the differences from the flat space, spin sum always leads to trace of some 
gamma matrices, the Fourier transformation of the fermion held is 

V’(x.i) = /'d 3 p^[e‘ px a p ,,X p , s (i)+e-‘ px 4, > »'p..(«)] 

J S 

= a“3(t) f d 3 p^[e ipx a P)S u P}S {t) + e~ ip ^/3l s v P}S (t)] . (Al) 

J s 

Namely, 

Xp, s (t) = a~^(t)u p>s (t), W P)S (t ) = a~^(t)v P)S (t). (A2) 


We define 


u s , P (r) 


I u +)P (t)^ s 
\U-, p {t)iI) s 




I V +tP (t)i/J-s 

\V-, P 


(A3) 


where ip + = \ and = 


0 

1 


stand for spin up and down, respectively, which make 


M SjP (r) and u SjP (t) connected by charge conjugation as v StP (r) = Cu^ (t) = A(td 7 °) r = 
C' 7 °ut I = C^°u*. Where 


C = 27°7 2 = 

is charge conjugation matrix. Then 


0 ia 2 
ia 2 0 


(A4) 


/ v+, P ( t)^- s 
\ v-, P (r)'tp-s 



I —icr 2 u*L p (t)iJj s 


(AS) 
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Therefore 


V +A T ) = V +A T ) X X 2x 2 = ^v+,p(r)^_ s ^ 


T 

-sr-s 


«-,pW = WpO) X 3-2x2 = J2 V ~’p( 



= -ia 2 u*_ (r)oi, (A6) 



= ^2«+, p M [ ” q | =io- 2 u+ iP (T)o- 1 . 
That is we have u± iP (r) = =Fi<T2'U^ p (T)<7i. We can now calculate the spin sum. 

s s' 

1 


a 3 (r')a 3 (T" 


^7ri[ u P 'A^) u p'A T ")]bcl2i[vpA^v p,®('r , )]da 


a 3 (r')a 3 (r") 1 

s,s' 


The trace is 


E tr[7 m u p / iS / (r')M p / iS /(r")7 n W,s( r ") v p,s (E 

s,s' 

ytr[7-n^ s ^0y,,(r'07 O 7%,(r / 0< s (-')7 O ] 
0 1 7 m \ I n +iP /(r')'0. 


S,S‘ 


E*[ 


CT m 0 J y M_ iP /(r')V'. 

0 cr n \ ( v +iP (t")iP- 




Cy(r'),f s ^ iP (T'))] 


cr n 0 / \ v- tP {r")ij)- a 

0 a m W M +iP /(r / )^''0jw+,p'(^") «+ >P '(' r/ )^ s ''0^-,p'(7' // ) 
a m 0 J \ v n_ i p/(r , )V ; 3''0^y,p/(r") n_ ! p^(r , )V' s ''0j^,p'(^") 

0 CTn \ f W,p('r")V , -s'0-W+,p(' r ') W,p(^ // )^- s '0-W-,p( r/ ) 


E tr t 


a n 0 / \ W-,p(r // )V’-s'0-W+,p(' r/ ) ^-.P^'OV’-sEEpO'') 


T 


while 


E^2 


S==b 




, 1 \ / 1 0 

0 , 1 + 1,0 = 

i/ Vo/ v o i 


(AT) 


Etr^p'.'Mv^OA.M^^')], (A8) 


(A9) 


(A10) 
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Therefore we have 


\ # /^p,^ 


tr['y m u p i tS > (r')upi :S ' (r")^ n Vp tS 

s,s' 

+, P ' (t')u ^ p , (r") M +iP / (t / )m 1 iP , (V 

-,P' ( r >+, P ' (t") w-,p' Cr>-,p' 0' 




\ v T,F v ^ N T' v 

v + ,p(r")vi,p(r') v +i p(r ,/ )vl,p(r') \ 

, x y-,p(r")4, P ( T/ ) v -,p( r ") v -, P ( r ') / 

tr / a m n_,p/(T>^p,(T") ^.^(rOM^p^r" 

y cr TO M +iP /(T')w+, p '('r") ffmM+ 1 p'(r> t _ iP ,(r"^ 

/ ffnU-.p^'Ou^p^) cr n u_ ) p(r ,, )ul i p(r / ) V 
\ ^n^Cr'Ov+.pCr') p (r') y 

tr[(T m ^^_,p/(r , )^^ 1 fI p/(r /, )o■ n ^;_ l p(r ,/ )^’+ I p(r , ) + u m u,_ 

+o■ m ^^+,p'( r, ) ^^1 f,p'( r,/ ) c^ r^W-,p(r")vlp(r , ) + a m M +) p/(r , )Ml i p,(r")crnW+ ! p v ,--p, „ 

tr[-a m M_ ! p/(r')M 1 fi p,(r")cr ri cr2M+ i p(r ,, )M^ ! p(rVl + crm“-,p'( r/ ) M - !P .(r ,/ )o'nO'2M-,p(r ,/ )M^ ! p 
+a m M + ,p/(r')M 1 fi p,(r ,, )cr n a2My i p(r")M+ ! p(rV2 - ^mU+, P '(r'W_ :pl (T ,, )a n a 2 u*_ p (T , ')ul tP (T , )al], 

(All) 


(JmU— : pi (t ) 


/(T / )Ml ip ,(r")cT ri 'i; +! p(T / >^p(T / ) 
l (r , >l : p(r / )] 


»1 


where we have used u± iP (r) = p (t)cti. The solutions are 


, V ( U V,p( T ) X i>s \ 

UsA T ) = . 

(A 12) 

y (p- (T ) u iiA r ) x i / 


with 


^+,p' ) l2x2'W/i,p(' r ); 

(A 13) 

W+.p'CO = (P- ^)%W, 

(A 14) 

where /i = 1 — ir, p, = 1 + ir and 


^ ) -2(2^/2 e2 ^ 1) (")’ 

(A 15) 

% ( X ) - 2(2^/2 e 

(A16) 

Hankel function goes as 



(A17) 
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in small x limit. Therefore in this situation, the mode functions behave like 


u p (x) ~ 
u^x) ~ 


T(/i)e 


(2vr) 2 
T(/i)e" 2 ¥ : (x\~ ir 
(2tt) 2 


(A18) 

(A19) 


Then we have 


s,s' 

= tr[—<T m w_ 


) u p', 

s yr"A n v P ,s{r")v p . 

t / 

//\ „ * / 


'+,p'( 

r )a n a 2 M +iP ( 

r )u 

( II 


\ T 

A T 

)a n a 2 u +p [T 

) u +,p 


) ff 2 


1 51 ^ ' It' ' ' “ Jl-' ' ' 5 JK ' ' ' 

■ / )w+ )P /(r")o-nO-2M; iP (T , >^ p (r , )(j| -(T m W +) p/(r> t _ !p ,(T // )o-n^2'ul, p (T")«+, p (T')cT2 
= tr[-a m (p' • a)u^ iP ,(r , )M; iP /(r")a n a 2 «*, p (r")(p ■ a) T u^ p {r')al 

+a m (p ■ a)u Ptpl (T')(p' ■ a) ] ul p ,(r")(T n a2{p • ^)*^, p (r")(p ■ a) T u^ p (r' 

+^m'W M ,p'(T / )M* )P /(T , > rl (T2'U* >p (r / > MjP (T / )o-| 

• ^) t «p, P '('T // ) CT n CT 2(p ‘ °Yu*^ p {t")u^ p {t')(t\} 

= tr [-(Tm{p' ■ a)a n a 2 {p • 

+a m (p' ■ a)(p' ■ <jy<J n a 2 (p ■ a)*{p • cr)^cr|^ ?p ,('r , )^ ;P '(' 7 " ,, )^p, P ( 7 ' ,, )'^/i, P (' 7_/ 

+ff m (7 n (72ff2 M M.'( r, ) M ftp'( r >ftp( T >M,p( T ') 

-a m (p' • a)V n a 2 (p • o') 1 

= tr[cr m (p' • a)a n (p • a)M A , p /(r / )M* iP /(r")M* iP (r")M A ,p(r / ) 

+ffmCT n Mp IP /(r , )^ IP /(r")M^ lP (r")Mp,p(r , ) 

+Om^ /Jl p'(T / )<, p '(T>; !P (r> ftP (T') 

(p' • °y°n(p ■ ^)Mp,p'(r , )M^ p /(r")uf )P (r")Mp,p(r')], 


+A 

where we have use 
(p-o-)(p-(y) ] = 


(A20) 


P3 Pi ~ IP2 
Pi + ip 2 P3 


P3 Pi - W 2 
Pi + W 2 -P3 


— (Pi +P2 + ^3)12x2, (A21) 


(P ' <7)*(P ' V) 1 = [{P' ■ <?)(?' ' ^) f ]* = 1 2x2, 


<72(p • o-) T a| = cr 2 (p • a)*al = -p ■ a, 


* A _ 


(A22) 

(A23) 
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while 


WM,pK r 0<, P '(O<p(^ / K,p(y) 

r(/x)r*(/i)r*(/i)r(/x) (- P w ,x ~ ir / —" N_ir/ — /N_ir 

(2^)8 

r 4 (/x) /- P v x " 2ir 


-pr 


—p r 

2 y V^T 

//\ —2ir 


—pT 


—pr 


(2tt) 8 V 2 / V 2 

u p,p'( 

WW'WfW /-A 

(2ir) 

r 4 (p) / -pV -pr” 

W)‘ 


/ / \ ir / / ft \ if / // \ zr / / \ zr 

• ~ x -p r \ / —pr \ / —pr x 






(A24) 


(A25) 


and 


Mp,p'(r , )n^ p ,(r")n^ p (r ,, )^ > p(r , ) 


r W r *M r * M r W ; - 2lf r 


(2tt) 

-e _27n r' r” 


/ /\ — zr 

-p r x 


|r(js)l 4 . „ — 2nrf—2ir//2ir 


(2t r) 8 

?W' (r'X,?/ (^>J» ftP ( T ') 


/ // \ ir 
-pT N 



-pr' 


«r 



(2tt) 8 

Finally, we get 


rMr-MrWM 27rr f-p'r'Y f-p'r"Y ir (-pr- ir 

(2tt) 8 6 

|r(/i )| 4 ^2nr ^_/2ir ^_f/—2ir 



-pr 


/ \ —zr 


—pr 


/ \ zr 


(A26) 


(A27) 


y^tr[7 m n p / ig /(r / )M p / iS /(r // )7 n n Pig (r // )n Pig (r / )] 

s,s' 

tr[a m (p' • cr)cx n (p • fT)Mp,p'(r')n piP ,(r'>* iP (r")n /i! p(r') 

+^m^nMp,p'(r')nt iP ,(r")n* p (r")wp i p(r / ) 

+^m^nMp,p'(r / )^ iP ,(r , >; p (r ,/ )np iP (r') 

+a m (p' • a)V n (p • a)7/p ! p/(r , )^ p ,(r ,, )Mt p (r")wp,p(r')] 
tr [a m (p' ■ a)cr n (p ■ 


(2tt) 


/ / \ —2ir 

-pr x 


+tr [a m cr n ]( 


|r(/i )| 4 

—27rr /—2zr _^//2z r 

(2vr) 8 


2 

^//zzr _|_ 


-pr 


—2ir 


|r(h)l 4 

(2t r) 8 


r 4 (/i) 

(2vr) 8 


-pV' 





-pr 


2ir 


^2nr ^/2ir ^f/—2ir\ 


(A28) 
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